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O ■ Abstract 

We show how to renormalize ^-derivable approximations in a theory with a fermionic field 
coupled to a self-interacting scalar field through a Yukawa interaction. The nonperturbative renor- 
malization concerns the self-interaction coupling of the scalar field which is renormalized via a set 
of nested Bethe-Salpeter equations for the scalar and fermionic four-point functions. We use this 
Q\ , information to construct explicit finite equations of motion in the symmetric phase. We work in 

the context of equilibrium quantum field theory and show that the renormalization can be carried 
out without introducing temperature dependent counterterms. 
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I. INTRODUCTION 



Formulation of quantum field theory in terms of the two-particle- irreducible (2PI) effective 
action has recently regained interest in the context of QCD. It leads for example to a good 
description of the entropy of the quark-gluon plasma from very high temperatures down 
to temperatures of the order of three times the critical temperature The 2PI 

effective action was introduced in the context of the non-relativistic many body problem 
and was then extended to quantum field theory in |7j]. It aims at parametrising 
physical quantities in terms of one- and two-point functions. For quantities dominated 
by fluctuations encoded in one- and two-point functions, an expansion of the 2PI effective 



action is then expected to show better convergence properties than a pertur 
expansion. This can be explicitly checked in equilibrium within scalar theories 
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that the same phenomenon could cure the poor convergence of 
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for particular thermodynamic quantities in QCD. 
Many other approaches aim at improving the convergence of perturbation theory using 
propagator resummation. 



theory for scalar fields 
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his is for example the case of the so-called screened perturbation 
Q 0r its generalization to QCD U H Q- These other 



approaches are based on a resummation of the leading contribution to the self-energy for 



the soft modes (Hard Thermal Loops [19|,|20j). However they also include these corrections at 
hard scales leading to spurious contributions, in particular artificial temperature dependent 
ultraviolet divergences. In contrast, the 2PI effective action treats soft and hard scales 
on a different footing thus preventing the presence of unwanted temperature dependent 
singularities. Finally the use of the 2PI effective action goes now beyond the description of 
systems in equilibrium since it has been successf ully applie d to describe late time dynamics 



of quantum fields out-of-equilibrium 
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30] and may then be 



useful in understanding equilibration properties of the system created in ultra-relativistic 
heavy-ion collisions. 

In order to do calculations within the 2PI framework, one has to devise approximations 
of the 2PI effective action. A particular class of approximations is obtained by considering 
diagrammatic truncations of the skeleton expansion [?| for the 2PI effective action. These 
approximations, known as ^-derivable 0, are important since they respect global symme- 
tries jj, 0, 0] and thus insure the presence of conserved quantities (for instance, energy 
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conservation is crucial in order to properly simulate the time evolution of quantum fields). 
There are however two major drawbacks when constructing ^-derivable approximations. 
The first concerns any renormalizable theory and has to do with the possibility or not to 
remove UV divergences within a given ^-derivable approximation. The second concerns 
gauge theories and how to preserve gauge invariance. The difficulties are related to the 
fact that a given truncation of the effective action only resums particular topologies and 
it is not clear how fundamental properties such as finiteness or gauge invariance manifest 
themselves within such a subclass of diagrams. Recent progresses have been made concern- 
ing renormalization in the case of scalar theories 

0, fl H Q H H Q Q, Q- 

For 

such theories, divergent topologies can be self-consistently renormalized by means of zero 
temperature counterterms. The question remains however totally opened to know what can 
be said about renormalization in the case of gauge theories. Two new features need to be 
understood: How to renormalize ^-derivable approximations coupling different type of fields 
and the interplay between renormalization and symmetry. We choose here to address the 
simpler one: As a first step we consider a situation with the same power counting rules as 
QED but without the difficulties related to the gauge symmetry. The aim is to see if the 
divergent topologies can be again self-consistently renormalized in a situation with two kind 
of fields. We thus consider the case of a fermionic field coupled to a self-interacting scalar 
field via a Yukawa interaction. We work in the context of equilibrium quantum field theory 
and discuss renormalization of the equations of motion in the symmetric phase. The broken 
phase will be dealt with in a subsequent work. 

Section II introduces the 2PI effective action together with the relevant equations involvec 



in the renormalization process. The result is indeed an extension of what is done in 
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371 ] . In particular, renormalization of singularities involving four scalar fields relies 



on the same Bethe-Salpeter equation. The new feature is that, the kernel of the latter 
arises itself from a second Bethe-Salpeter equation for a four-point function involving only 
fermionic fields. This set of nested Bethe-Salpeter equations is the key to disentangle the 
coupling between scalar and fermionic degrees of freedom in the analysis of UV divergences. 
In section III, we use this information to renormalize the equations of motion in the one-loop 
approximation. Finally, Section IV is devoted to the generalization to any loop order. We 
provide a systematic rule for renormalization. We conclude by saying a few words on the 
difficulties met when extending these results to Abelian gauge theories. 



3 



II. DIVERGENT TOPOLOGIES OF THE 2PI EFFECTIVE ACTION 



We consider a fermionic field ip coupled to a self-interacting scalar field tp (if 4 theory) via 
a Yukawa interaction, in four dimensions. We consider for the moment the unrenormalized 
(bare) action which reads 



The time integral runs from to —if3 where (3 = 1/T denotes the inverse temperature. 
This theory is perturbatively renormalizable: The couplings in the action are dimensionless 
or, equivalently, the superficial degree of divergence 5 of any Feynman diagram does not 
depend on the number of vertices but only on the number of scalar (E v ) and fermionic (E^) 
external legs: 5 = 4 — E v — (3/2)E^. The theory has to be regularized so that one can 
write meaningful expressions in auxiliary steps before renormalization. In this paper, we 
use dimensional regularization which amounts to working in dimension d = 4 — 2e. In that 
case, the couplings acquire a dimension which we parametrize by a renormalization mass 
scale fi: g —> fi e g and A — ► /x 2e A. 

We use the Minkowski metric (+, — , — , — ) and work in Fourier space. We denote 4- 
momenta by capital letters Q = (qo,q). At finite temperature in imaginary time, the fre- 
quency integrals or sums are performed over the imaginary axis Q = (qo = iu),q). We use 
the following notations: 



A. 2PI effective action and equations of motion 

In the symmetric phase, the mean values of both fields vanish and in order to study 
thermodynamic properties of the system, it is enough to consider the 2PI effective action as 
a functional of the full propagators S and D. It is given by the expression 




(1) 






-i 



-Tr ILD + $[S, D; g, A] . (3) 
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The trace symbol Tr includes (d — l)-momentum integrals, Matsubara sums, and sums 
over spin indices. The self-energies £ and n are related to the propagators by the usual 
Schwinger-Dyson equations: 

S~l{P) = -f + m + E Q/3 (P) , D-\K) = K 2 + M 2 + U(K) . (4) 

The functional $[S, D; g, A] is given (up to a minus sign) by the sum of zero-leg two-particle- 
irreducible diagrams drawn with propagators S and D, and vertices g and A. In Fig. ^ we 
represent the first diagrams contributing to $. The functionals $ and r 2PI depend a priori 
on the regulator e and on the renormalization scale /i even if we do not make it explicit in 
our notation. 




FIG. 1: First contributions to <&[S, D; g, A]. Dashed and solid lines denote scalar and fermionic 
propagators respectively. The vertices are bare, tree level vertices. 



The physical propagators are given by the stationary condition on r 2PI [S, D] which is 
suitably translated into a set of equations of motion: 

These equations state that, in order to obtain the stationary self-energies, one has simply to 
remove respectively a fermionic or a scalar line from the diagrams in $. By evaluating the 
2PI effective action @3J) at its stationary point, one obtains an expression for the pressure of 
the system in terms of the full propagators S and D. Other thermodynamic quantities can 
be obtained by taking derivatives of the pressure. 

^-derivable approximations are defined by selecting a class of diagrams contributing to $ 
in Eq. Q and changing the equations of motion (j3J) accordingly. If the diagramatic trunca- 
tion of $ is done according to the number of loops we refer to it as the 2PI-loop expansion. 
Further approximations can be build on top of the 2PI-loop expansion such as the ones used 
in jl, Q, [| or in j3, 13, Q (see also J2 Q])- They have the advantage to be easier to 
implement numerically than a full self-consistent calculation. However the first approach is 
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not systematic and the second is based on an expansion which leads to unbalanced tempera- 
ture dependent singularities. In contrast the 2PI-loop expansion is systematic and as we will 
show can be renormalized without introducing artificial counterterms. Before dealing with 
general 2PI-loop app roximations in section I1VI we illustrate the renormalization procedure 
58| for which the equations of motion read 



at one-loop level 



rK 1 ) 

S Q/9 (P) = -g 2 / D(Q)S a p(Q + P), 

JQ 

U(K) = -A / D(Q) + g 2 / tx[S(Q)S(Q + K)j . (6) 

J Q J Q 

The notation includes finite temperature contributions via Matsubara sums. The trace 
tr is taken over spin indices and we have taken into account the minus sign arising from the 
fermion loop. These equations are represented diagrammatically in Fig. El 

/'~^ 

£ =- > 'S > n = - 'C^y ' ' 

FIG. 2: Equations of motion at one-loop level (the lines represent full propagators). 



B. Iterations and divergences 

The aim of this section is to reveal some of the topologies which play a role in the process 
of renormalization at one-loop level. This simple example is indeed very fruitful since the 
structure of the equations that it leads to plays a general role (see section HVjl . Most of the 
ideas are those already described in (stJ ■ The strategy to study divergences in the equations 
of motion is to iterate them, to express formally their solution in terms of an infinite sum of 
perturbative diagrams and to ask whether it is possible to tune the counterterms in order 



to absorb the singularities |59|. In order to do this consistently, one should proceed order 
by order in A and g. Here we just give examples of the divergent topologies that one reveals 
by this procedure and claim that there are no other ones. Later we prove this statement by 
providing renormalized equations of motion. 

At one-loop level we will be only concerned by singularities involving four scalar fields 
(contributing to the scalar coupling counterterm: SX). Field strength singularities are dealt 

6 



with when treating the asymptotic equations of motion (see section UlI C|) . Mass countert- 
erms do not play a role since we restrict ourselves to the massless theory in dimensional 
regularization j60( • The Yukawa coupling counterterm 5g is treated in section II VI when 
dealing with higher loops where its renormalization becomes important. At the level of the 



one- loop approximation, Sg = 0. Things are different for the sea 



ar coupling counterterm S A 



which receives contributions to all orders. As was shown in [3jJ in the case of the pure ip 4 
scalar theory, an analysis of divergences based on BPH (Bogoliubov-Parasiuk-Hepp) boxes 
I52I [53I . with full propagators misses a large amount of coupling singularities. The idea 
is that if one views the propagator as build from a resummation of diagrams, one creates 
new boxes which are not taken into account by the BPH procedure applied to resummed 
diagrams. In Fig. El we give an illustration of such a phenomenon in the scalar ip A theory. 
We will show that these singularities can be accounted for by a (temperature independent) 
redefinition of the scalar coupling counterterm in front of the tadpole diagram (see for in- 
stance Fig. |2I). As we will show, this redefinition is governed by a Bethe-Salpeter equation 
0. 





(i) 

(2) 




(2) 





(1) 




FIG. 3: The BPH procedure on resummed diagrams misses singularities. This is simply because 
the two operations drawing boxes (1) and iterating diagrams (2) do not commute as illustrated 
in the figure. The first line shows the coupling singularities which are accounted for by a BPH 
analysis regardless of the content of the propagator. The second line shows what happens if we 
first look inside the propagator and then apply the BPH analysis. There are clearly more divergent 
topologies in the second approach. 



In our present example, the singularities appear in the equation of motion for the scalar 
propagator after iterating the coupled equations of motion a certain number of times. By 
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iterations we mean that, knowing the solutions of these equations X*™) and Il( n ) to a given 
accuracy, we can generate new contributions by plugging the related propagators = 
Sq 1 + £ (n ) and = G 1 + into the equations of motion and expanding perturbatively 
according to 

s in) = s -s e^s + s e^s e^s + ... 

G {n) = Go - GqU^Gq + G U^G U^G + ... (7) 

The simplest such a divergence appears by considering the Oth iteration of the equation of 
motion for the fermionic propagator (i.e. the perturbative one-loop contribution) and plug- 
ging it into the equation of motion for the scalar propagator. After opening the scalar line 
(opening a line corresponds, in the BPH language, to drawing a box which pulls this line 
apart) one ends up with the diagram in the right hand side of Fig.HJ which is logarithmically 
divergent according to its superficial degree of divergence (5 = 0). A simple generalization 



~ •--<>- ' 



FIG. 4: UV singularity generated by plugging the fermion self-energy to lower order into the 
equation of motion for the scalar propagator. The lines represent perturbative propagators. 



of this construction consists in first iterating the equation of motion for the fermionic prop- 
agator by keeping the scalar propagator fixed to its free value So- We generate the rainbow 
diagrams of Fig.|5]which, once plugged into the equation of motion for the scalar propagator, 
lead to the singularities of Fig. El which generalize the singularity in Fig. HJ The first two 
steps are the same than in Fig.HJ After iterating (step 1), one selects a scalar line to isolate 
a singularity with four scalar legs (step2). Step 3 only resums the free fermionic propagators 
Sq into the full propagator S. The final result of this procedure is the general structure on 
the right of Fig. H which generalizes the one in Fig. 0] by summing all the possible ladders 
made of scalars rungs. 

One can finally repeat the same construction by considering corrections to the scalar 
propagator. If we consider for example the scalar self-energy correction that we discussed in 



S 



FIG. 5: Rainbow diagrams generated by iterating the fermion equation of motion with fixed scalar 
propagator (S = So). 
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(cl) 



(dl) 
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r*~t 












FIG. 6: The rainbow diagrams modify the singularity of Fig. All fermionic lines represent 
perturbative propagators So- Only in the last step (3) these are resummed into S. 



the previous paragraph and take it into account in the generation of rainbows, one obtains 
the divergences depicted in Fig. which are nothing but a concatenation of those depicted 
in Fig. H3 The general divergent topology contributing to 5X is then the one depicted in the 
first diagram of Fig. |HJ The tadpole diagram brings trivial modifications to this topology 
(see the second diagram of Fig. [HJ). We claim here that these are the only singularities 
generated by the equations of motion which contribute to 5X. We give an algebraic proof of 
this statement when renormalizing the equations of motion in section HTT1 
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FIG. 7: Structures such as the ones depicted in Fig. are concatenated through the iteration of 
the equations of motion. 




FIG. 8: General structure of scalar coupling singularities in the one-loop approximation. 



C. Kernels and Bethe-Salpeter equations 

Let us for the moment analyze in further details the topologies described in the previous 
section. We show that it is possible to simply encode them using nested Bethe-Salpeter 
equations. To this purpose, it is convenient to introduce the following 2PI kernels (the 
particular ordering of labels is introduced for later convenience): 



^ SSpaSSys ' ^ SSpaSD ' 5D5S aP ' w SD6D ' 

(8) 

The numerical factors are introduced so that the kernels represent (up to a minus sign) the 
sum of related diagrams. In Fig. El we represent them at one-loop level. In order to generate 
the singularities of Fig. |H1 we start from the first kernel which in the one-loop case reads 
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FIG. 9: Basic kernels in the one-loop approximation. 



We iterate it via a first Bethe-Salpeter equation for the four-point function with four 
fermionic legs: 

T^(P,K) = A^(P,K)- [ A^(P,Q)M(Q)T^(Q,K), 

Jq 

= A^(P,K)- [ T^(P,Q)M(Q)A^(Q,K). (10) 
Jq 

with M( Q , i g) ) ( 7 5)(<5) = S a ~ f (Q)Ss/3(Q)- The indices in (0) have been defined so as to be able to 
use a generalized matrix product in (|10|) (see Appendix El . This Bethe-Salpeter equation 
resums all the contributions of the type given in Fig. [TQJ We now use the kernels A^ v and 




FIG. 10: First Bethe-Salpeter equation. 



A w together with T^, to define a new kernel as follows: 

A W (P,K) = A VV (P,K)+ [ A^(P,Q)M(Q)A^(Q,K) 

Jq 

-If A^(P,Q)M(Q)T^(Q,R)M(R)A^(R,K) , 
Jq Jr 

= A + V / S\P + Q)M(Q)S(Q + K) 
Jq 

- 4# 4 / [ S\P + Q)M(Q)T^(Q,R)M(R)S(R + K), (11) 
Jq Jr 

where we used the one-loop expressions A W (P, K) = A, A^(P, K) = —2g 2 S(P + K) and 
= A\ . The upperscript t denotes transposition in the generalized index space (see 
Appendix A). Note that we do not need to use the trace symbol since this is accounted for 
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FIG. 11: New kernel. The factors of 4 present in equation (|11|) and the signs arising from fermion 
loops are not represented since they are included in the diagrammatic representation. 



by the generalized matrix product. This kernel, which we represent in Fig. ^2 generates 
all the coupling singularities that we discussed in the previous section via a second Bethe- 
Salpeter equation for the four-point function with four scalar legs: 



T W {P,K) = A w (P,JO-- / A iptp (P,Q)D 2 (Q)T iptp (Q,K) 

J Q 

= A w (P,K)-± I T (p(p {P,Q)D 2 (Q)A^(Q,K) 



(12) 



This Bethe-Salpeter equation [62] represented in Fig. contains UV di verg ences and is 



renormalized in the very same way than in the case of scalar theories 
analysis is given in section ITll Bl 



33, |37|1 . A detailed 
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FIG. 12: Second Bethe-Salpeter equation. It encodes the coupling singularities of the equation of 
motion for the scalar propagator. 



III. RENORMALIZATION 

We now prove that the scalar coupling singularities present in the equations of motion 
are indeed those encoded in the Bethe-Salpeter equation (fI2|). To this aim we show that, 
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once the equations of motion are expressed in terms of the finite function T W (P, K), there 
are no other divergences. In this section, we illustrate this point at one-loop level. Higher 
loops are discussed in section HVI 

A. Renormalized quantities and counterterms 

In order to perform renormalization, it is convenient to express the equations of motion in 
terms of renormalized (finite) parameters and counterterms (infinite). To make an explicit 
difference between renormalized and unrenormalized quantities, we label the latter by "b" 
(bare). We then rescale the fields according to ip b = Z^ 2 ip and tpt, = Z^ 2 tp so that the 
action becomes 

S = j d d x \z^(i<f) - m b )^ - g h Z^Z l J 2 i,^ + ±Z v {d<p) 2 - \m^ 2 - ^A b Z> 4 } . 

(13) 

Next we introduce renormalized parameters and counterterms by setting 

Z f m b = m + 6m, Z^M 2 = M 2 + 5M 2 , g h Z^ 2 = g + 6g , X b Z 2 = A + <5A , (14) 
and write the action in the following way: 

S = J d d x j^(z0 - m)V - gW>H> + \(dv? - - \^ 

+i5Z^dip - Sm^ip - 6g^<p + ^6Z v (d<p) 2 - ^5M 2 <p 2 - ^A^ 4 J , (15) 

where we have defined 5Z V = Z^ — 1 and 8Z^ = Z$ — 1. We also define the renormalized 
propagators S = ZT S b and D = Z~ l D b together with the renormalized self-energies (we 
will indeed show that these are the objects which can be made finite) S(P) = S~ 1 (P)+]f , — m 
and U(K) = D~ X (K) - K 2 - M 2 . We thus have 

E(P) = -f>6Z^ + 5m + Z^ b (P), U(K) = K 2 5Z^ + 5M 2 + Z v ,U h (K). (16) 

It is then easy to express Z^Y, b and Z^\. b in terms of renormalized quantities: In the 

expressions for S b and lib, one does the replacements 5b — > S, D b — > D, g b — > g + 5g and 

Ab — ► A + SX. Thus in the one-loop approximation, the equations of motion read 

r(T) 

E(P) = -(g + Sgf D(Q)S(Q + P)-fSZ^ + 8m, (17) 

J Q 

1 r{T) r{T) 

U(K) = - (A + 5X) / D(Q) + (g + Sgf / S\Q)S(Q + K) + K 2 5Z^ + 5M 2 . 

2 JQ JQ 
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These equations are UV divergent. We shall show that it is possible to absorb consistently 
all the UV singularities in the counterterms present in (|17|). In what follows, and to simplify 
the discussion we study the massless case both at zero and finite temperature. Since we 
use dimensional regularization 5m = and 5M 2 = 0. Other regularizations, such as using 
an explicit cut-off, generate linear and quadratic divergences which can be simply absorbed 
in mass counterterms. As already mentioned, the Yukawa coupling is not renormalized at 
one-loop level. We then set 8g = 0. 

Finally the renormalized kernels are obtained by factoring out the appropriate Z factors: 
= Z^A^, A^ = Z^Z^A^ and A w = Z 2 A h w . They are obtained from the bare kernels 
by doing the replacements — > S , D h — > D , g h — > g + 5g and Ab — > X + 5X. Note that these 
kernels do not need to be finite since, in the exact theory, they do not correspond to any 
physical four-point function. It turns out that A^ and A^ can be made finite but not A w 
since it contains the counterterm 5X which is used to renormalize the four-point function 
T w = Z^r^ . At one-loop level, we have 

A^ = -g 2 D , = -2g 2 S , A w = A + 5X . (18) 

The equations defining = Z^T^, A w = Z 2 A h w and T w = Z 2 Y h w are exactly those 
given in (|1U|). (jlljl and (|12p after one replaces bare by renormalized quantities. 



B. Nested Bethe-Salpeter equations 

The first step in the renormalization procedure is to renormalize the set of nested Bethe- 
Salpeter equations. The Bethe-Salpeter equation (jlUJI defining is finite. In the one- 
loop case, this follows by simple inspection of the diagrams which are resummed by the 
equation. The only subdiagrams involving loops have four fermionic legs and thus have a 
negative superficial degree of divergence (8 = —2). When used to build the kernel A w , 
the function generates logarithmic singularities since A w has four scalar legs (5 = 0). 
However the singularities are all overall singularities since there is no divergent subdiagram 
(all subdiagrams have at least two fermionic legs and thus yield a negative superficial degree 
of divergence 5 < —1). 

The Bethe-Salpeter equation defining T w is the same than the one described in [37]]. Its 
renormalization thus follows the very same lines. It amounts to adjusting the scalar coupling 
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counterterm 5X [63] present in A W (P, K) and that originates from the tadpole diagram. The 
finite equation for T W (P, K) reads: 

- i y {a w (p, q) - a w (p„ Q) }p 2 (g)r w (g, #) 

- ^r w (p„g)p 2 (g){A w (g ! ^)-A w (g,^)} ! (19) 

where (P^fC) denotes the renormalization point. The finiteness of this equation follows 
from the fact that A W (P, g) — A W (P*,Q) is finite (there is only an overall logarithmic 
divergence in A w ) and the property A W (P, Q) — A W (P*, g) ~ 1/g for large Q and fixed P 
and P*. In pure y? 4 theory, this property is related to the 2PI structure of the kernel A w (37 1. 
Here the kernel A w is not fully 2PI since it arises from a Bethe-Salpeter equation which 
is two-particle-reducible (2PR) with respect to fermionic lines (see Fig. ITT)j) . However this 
does not prevent the asymptotic property of A w to remain true as we discuss in Appendix 

m 

C. Asymptotic self-energies 

An helpful quantity in the discussion of UV singularities is the dominant asymptotic 
piece of each of the self-energies (37 1. We denote them respectively by Si(P) and n 2 (A'). 
The labels are used to recall that, for large momenta Si(P) ~ P and I\.2{K) ~ K 2 up to 
logarithms. These asymptotic pieces satisfy simplified equations of motion obtained after 
removing all the mass scales (but the renormalization scale fi) in the initial equations of 
motion (fTTjl . For the one- loop case, we obtain 

Sx(P) = -g 2 [ D 2 (Q)S 1 (Q + P)-f5Z^, 
Jq 

U 2 (K) = l(X + SX) [ D 2 (Q)+g 2 [ S t 1 (Q)S 1 (Q + K) + K 2 5Z V , (20) 
1 Jq Jq 

with the asymptotic propagators S'i(P) = l/(-f + Si(P)) and D 2 (Q) = l/(Q 2 + n 2 (Q)). 
Note that, since we are working in the massless case, the asymptotic equations are nothing 
but the equations of motion in the vacuum. 

The renormalization of the asymptotic equations of motion is done by means of the field 
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strength counterterms 5Z^ and 8Z^. The choice 

SZ 4 , = -g 2 -^J^D 2 (Q)S 1 (Q + P) 
6Z V = -g*^Lf Sl(Q)S 1 (Q + K) 

J (J 



p* 

(21) 



K t 

leads to finite equations [64]. Note that the field strength counterterms only depend on the 
dominant asymptotic pieces of the self-energies. Consecuently they are those needed in the 
initial equations of motion (fTTj) . The value of the coupling counterterm 5\ does not play a 
role since in dimensional regularization J^D 2 (Q) = 0. 

Once the dominant pieces are known, we only need to focus on the subleading pieces 
£_i(P) = E(P) - Ei(P) and Tl (K) = U(K) - U 2 (K). The labels are again related to the 
large momentum behavior: XL^P) ~ 1/P and Uq(K) ~ 1, up to logarithms. For example, 
for the fermionic self-energy, the subleading behavior arises by fixing the momentum in one 
of the lines of the diagram in Fig. El The subgraph thus obtained is of degree of divergence 
5 < and according to Weinberg's theorem |55j decreases at least as 1/P. In order to obtain 
equations for S_i(P) and H (K), we use expansions around the asymptotic propagators: 

S(P) = S^P) + 5S(P) , D(K) = D 2 (K) + 6D(K) , (22) 

with: 

5S(P) = -SxiP^-^S^P) + S r (P) , 5D{K) = -D 2 (K)U (K)D 2 (K) + D r {K) , (23) 

and S r (P) ~ 1/P 5 and D r (K) ~ 1/K 6 for large momenta. In the massless case, at finite 
temperature, it is important to keep an infrared regulator which is provided by thermal 
corrections to the propagator. It is then more convenient to think of an expansion of Si and 
D 2 around S and D respectively: 

S 1 (P) = S(P)-5S(P), D 2 (K) = D(K) — 5D(K) , (24) 

with: 

5S(P) = -S(P)S_x(P)5(P) + S' r (P) , 8D(K) = -D(K)U (K)D(K) + D' r (K) , (25) 

with S' r (P) and D' r (K) keeping the same asymptotic behavior than S r (P) and D r (K) re- 
spectively. 
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D. Equations for E_i and IIo 

Before writing the equations for S_i and IIo, it is convenient to separate the explicit 
thermal dependence in the equations of motion. To this aim we perform the Matsubara 
sums. In the one-loop case, the thermal separation is easily achieved (see Appendix : 

S(P) = -g 2 [ D(Q)S(Q + P)-P5Z^ 
Jq 

- g 2 [ a^Q)S(Q + P)+g 2 [ D(P + Q)a^{Q) , 
Jq Jq 

U(K) = Ux + 5X) [ D(Q)+g 2 [ S t (Q)S(Q + K) + K 2 5Z^ 
2 Jq Jq 

+ \ (A + 8X) [ a v (Q) - 2g 2 [ a\{Q)S{Q + K) . (26) 
z Jq Jq 

The first lines of each equation are the same functionals of S and D that we would have 

written in the vacuum. They contain however an implicit thermal dependence via the 

propagators. The other terms (explicitly depending on temperature) are expressed in terms 

of (T^iQ) and (Tip(Q) which are defined as follows: 

M<2) = e(go)n|g |P^(g , q) , M<2) = ^o^mM*, q) , (27) 

with e(qo) denoting the sign function, ni ?0 i and A^| go | the fermionic and scalar thermal factors 
(they vanish as T — >• 0) and p^,(qo, q) and p^{qo, q) the fermionic and scalar spectral densities 
respectively (see Appendix EJ) . Finally Q denotes Minkowskian momentum. Note that if it 
were not for the counterterm 5X multiplying the explicitly temperature dependent tadpole 
integral, all the explicitly temperature dependent pieces would be finite. For higher loops 
one would need to exploit the techniques developed in [56] in order to separate explicit 
thermal dependences. 

We now derive equations for the subleading pieces E_i and n . We start from the 
asymptotic equation of motion for the fermion (f2ll|) and expand the propagators S\ and D 2 
by means of ()24|): 

E X (P) = -g 2 [ D(Q)S(Q + P)-f5Z i , 
Jq 

+ g 2 f 5D{Q)S(Q + P)+g 2 [ D(P + Q)5S(Q) -E r (P) . (28) 
Jq Jq 
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We subtract this equation from (|26|) and obtain an equation for E_i: 

= -<? 2 / 5p(Q)s(q + p)-s 2 / J D(p + Q)55(g) + s r (p) 

./Q JQ 

- g 2 f <r<p(Q)S(Q + P) + g 2 [ D{P + Q)a^Q) , (29) 
Jq Jq 

where E r (P) is a finite function decreasing as 1/P 3 as it can be checked on the explicit 
formula by using power counting and Weinberg's theorem: 

E r (P) = g 2 I 5D(Q)SS(Q + P) . (30) 
Jq 

In order to reveal the self-consistent structure of equation (J2HJ), we isolate the dominant 
contribution to 5S in (j23|) : 

E_i(P) = -g 2 [ 5D(Q)S(Q + P)-g 2 f a v (Q)S(Q + P) 
Jq Jq 

+ g 2 [ D(P + Q)M(Q)^ 1 (Q) + K(P)- (31) 
Jq 

E^,(P) is a finite function, still decreasing as 1/P 3 : 

E' r (P) = E r (P) - g 2 [ D(P + Q)S' r (Q) + g 2 f D(P + Q) a^Q) . (32) 
Jq Jq 

The last term seems to contradict the 1/P 3 behavior since the leading contribution at large 
P is —g 2 D(P) Jq cr^(Q) ~ 1/P 2 . However, as we show in Appendix O the integral of cr^ 
vanishes. The asymptotic behavior is thus ~ 1/P 3 . In Eq. ()31|) we have kept the dominant 
contributions to the asymptotic behavior 1/P of E_i. Note that these pieces are finite 
(this is expected since the equation of motion for the fermionic propagator contains only 
a field strength singularity which is entirely accounted for by the asymptotic equation of 
motion). However these dominant pieces, when plugged in the equation of motion for the 
scalar propagator, generate new singularities which need to be removed. The equation of 
motion (pTTj) is represented in Fig. O 

We now proceed in exactly the same way with the equation of motion for the scalar 
propagator and obtain: 

U (K) = l(X + 6X) f 6D(Q) + l{X + SX) f_a v (Q) 
1 Jq 1 Jq 

- 2g 2 [ S t (K + Q)M(Q)E^ 1 (Q)+W r (K), (33) 
Jq 
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FIG. 13: Equation of motion for S_i. The dotted lines in the first two diagrams represent the 
convolution of 5D(Q) and S(P + Q) or a v (Q) and S(P + Q). They do not represent propagators. 

which is represented in Fig. El n',(i^) is given by: 

U' r (K) = -g 2 [ 5S t (K+Q)5S(Q)+2g 2 [ S t (K+Q)S>(Q)-2g 2 f a\,(Q)S(K+Q) . (34) 
Jq Jq Jq 

Again by power counting and applying Weinberg's theorem we obtain that U' r (K) is a finite 




- + n r 



FIG. 14: Equation of motion for Ilo. The dotted lines in the first two diagram represent the 
convolution of 5D(Q) or cr^(Q) with 1. They do not represent a propagator. Symmetry factors 
and signs arising from fermion loops are included in the diagrams. 

function which goes like 1 / K 2 for large momentum. Unlike the fermionic gap equation, there 
are logarithmically divergent contributions in equation (J33)) . For the tadpole 5 = 4 — 4 = 
and for the fermionic insertion 5 = 4 — 3 — 1 = 0. These logarithmic divergences are those 
hidden in the Bethe-Salpeter equation for T w . The game to play is then to try to use the 
information on the four-point function T w to get rid of the singularities. 



E. Iterating the equation for E_i 

We now work on equations ()31|) and (|33j) or Fig. ^] and El We start manipulating the 
equation of motion for the fermionic propagator. We notice that the kernel appears 
in the third term of We can thus replace it by T^, by means of the Bethe-Salpeter 
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equation (JTUJ). This is diagrammatically shown in the first line of Fig. E3 We obtain 
S_x(P) = -g 2 [ 5D(Q)S(P + Q)-g 2 [ <j v {Q)S{Q + P) + E' r (P) (35) 

J Q JQ 

- [ T Hl (P,Q)M(Q)^ 1 (Q) - [ [ T^(P,K)M(K)A^(K,Q)M(Q)^ 1 (Q). 

JQ JQ J K 

We use the equation of motion (J31j) or the Fig. EH to recognize part of the integrand in (J35j) 




FIG. 15: Using the first Bethe-Salpeter equation and constructing part of the kernel A. 

(we hide irrelevant terms in E"): 

S_i(P) = -g 2 [ 5D(Q)S(P + Q)-g 2 f <x v (Q)S(Q + P) + £?(P) 

JQ JQ 

+ 9 2 [ [ T^(P : Q)M(Q)S(Q + K)8D(K) 

JQ JK 

+ g 2 [ fr^(P,Q)M(Q)S(Q + K)a tp (K), (36) 
Jq J k 

with: 

K(P) = K(P) - I r^(P, Q)M(Q)K(Q) ■ (37) 

JQ 

The properties of £' r '(-P) in the UV are those of H' r (P). This manipulation is also represented 
in Fig. ED Note that we have constructed part of the dressed kernel A w that we need to 
generate the UV singularities. 
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F. Iterating the equation for IIo 



Let us now consider the equation for the scalar and plug the expression for S_ithat 
we have just obtained. Again up to finite quantities which we hide in IT", we obtain 

no(ff) = \ (X + 5X) f 5D(Q) + ~(\ + 5\) fa v (Q)+W;(K) 

J Q J Q 

+ 2g A ! ! S\K + Q)M(Q)S(Q + P)5D(P) 
Jq Jp 

+ 2g A [ [ S\K + Q)M{Q)S{Q + P)a^P) 
Jq Jp 

- 2g 4 f f [ S t (K + Q)M(Q)T^(Q 1 K)M(K)S(K + P)5D(P) 

Jq Jk Jp 

- 2g 4 f [ f s t (K + Q)M(Q)T^(Q,K)M(K)S(K + P)a^P), (38) 

Jq Jk Jp 

with: 

U';(K) = U' r (K) - 2g 2 [ S^K + Q)M(Q)K(Q) , (39) 

Jq 

which again shares the same properties as U' r (K) in the UV. The manipulation is shown in 
Fig. ED Equation (|3*Hj) may look complicated but we have indeed reconstructed the kernel 
A w (see Eq. ITTj) and finally 

U (K) = l -j A tpv (K,Q)6D(Q) + ^ jk w {K,Q)a v {Q) + Ii l :{K). (40) 

J Q J Q 



G. Renormalization 

We recover the situation which we had with scalar theories [3j|. If we no" 
Bethe-Salpeter equation (fT^j). we get 

IL {K) = i ! T w (K,Q)5D(Q)+ l - I [v w (K,Q)a v (Q) +K(K) 

J Q J Q 



Q 

I 

'Qjp 



+ \ I I T W (K,P)D 2 (P)A W (P,Q)5D(Q) 




\\l T w (K,P)D 2 (P)~A w (P,Q)a v (Q) , (41) 

Q Jp 
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FIG. 16: Constructing the kernel A, 



The equation of motion (|4()jl allows then to write 

U (K) = l -l T w (K,Q)5D(Q)+ l - I [r w (K,Q)a v (Q] +K(K) 

J Q J Q 

+ i Jj w (K J Q)D 2 (Q)[u (Q) -KiQ)}. (42) 

Using the explicit expression for 5D(Q) given in (|25jl . we check that the potentially divergent 
pieces compensate leaving the equation equation: 

Uo{K) = K(K) + k [. r vp( K >QK(Q) + \ I r w (K,p){D' r (p) - D\p)n>;.(p)} . (43) 

Using the asymptotic behavior of each of the functions apperaring in the integrands, it is 
easy to check that this equation does not contain UV divergences. This ends the proof of 
renormalizability. 



IV. HIGHER LOOPS 

We now show that the general proof follows the same lines than that of the one-loop case. 
The main idea is that the equations of motion generate the nested Bethe-Salpeter equations 
which we need to account for the scalar coupling singularities in n . However we first need 
to describe the Yukawa coupling singularities in the general case. 
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A. General approach on the 2PI effective action 



To understand how the renormalization of the Yukawa coupling works, we first apply a 
BPH procedure on resummed diagrams. To this aim we make the counterterms explicit in 
the 2PI effective action. One can easily check that this amounts to replacing bare propagators 
by physical propagators, bare couplings by finite couplings plus counterterms and adding to 
$ two extra diagrams carrying the field strength counterterms. We show $ in Fig. El where 
we have omitted the renormalized couplings. The idea is now to apply a BPH procedure 

s g 

-* =:'"' '|j8Z + ^^{SZ^ + ( Y ) + 5g^^^5g + 8X|" + Sg^*^ ISg 

FIG. 17: $ functional expressed in terms of physical parameters and counterterms. 



on resummed diagrams. This means drawing boxes which contain divergent diagrams and 
asking whether there is a structure present in $ to account for this singularity. We have 
three kind of boxes depending on the number of external legs: 2-point boxes are related 
to field strength singularities, 3-point boxes to the Yukawa coupling renormalization and 
4-point boxes to the scalar coupling renormalization. 

Let us start drawing 2-point boxes as shown in Fig. El Because of the 2PI structure, 
these boxes contain all the lines of a given diagram but one. There are two diagrams 
waiting to absorb these singularities, namely the first two diagrams in Fig. El carrying the 
field strength counterterms. This procedure correctly accounts for 2-point singularities and 
even gives a formula for the field strength counterterms in terms of the propagators 
which generalizes Eq. (|21j) . 

We continue now with 3-point boxes as shown in Fig. El We see that we have to go to 
the last diagram to be able to draw such a 3-point divergent box. Because we are working 
withing ^-derivable loop approximations, we know that there is a topology present in $ to 
absorb this singularity, namely the fourth diagram in Fig. El This procedure fixes the value 
of the counterterm Sg in this diagram. Note that at this order of approximation 5g = for 
the last diagram in Fig. El This is nothing unusual. When renormalizing at a given order 
in perturbation theory we use different approximations of the same counterterms in order 
to construct a finite expression. In particular in the higher loop diagrams 5g = 0. From 
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C> O C> (^y 

FIG. 18: BPH procedure to renormalize 2-point singularities. The BPH procedure on resummed 
graphs is enough to account for this kind of UV singularities. 

this analysis it is clear that not all ^-derivable approximations are renormalizable. If one 
chooses a given diagram to appear in one has to include all the topologies related to it 
by the BPH procedure (^-derivable loop approximations fulfill this property). Furthermore 
in order to renormalize such an approximation, one has to allow the counterterms Sg to be 
different from one topology to another. This procedure accounts for the divergences that 
need to be absorbed in the Yukawa coupling renormalization. 




FIG. 19: BPH procedure to renormalize 3-point singularities. The BPH procedure on resummed 
graphs is enough to account for this kind of UV singularities. 

Let us finish by mentioning that in the case of 4-scalar singularities, the same remarks 
made at the level of the equations of motion still hold here. Namely, the BPH approach 
applied to diagrams with full propagators misses coupling singularities. However all these 
extra singularities are, as in the equations of motion, accounted for by a modification of the 
counterterm in front of the tadpole. This extra contribution accounts for all the coupling 
singularities generated by the Bethe-Salpeter equation (see Fig. EH). To account for all these 
singularities, we only need to build the correct kernel A. To that aim, in the next section, 
we look into the equations of motion at any loop order. 
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FIG. 20: BPH procedure to renormalize 4-point singularities. The BPH procedure on resummed 
graphs is enough to account for this kind of UV singularities. This procedures misses singularities. 
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FIG. 21: The first diagram represents the kind of scalar coupling singularities 5A BPH accounted 
for by a BPH procedure on resummed diagrams. This procedure misses singularities. However all 
these singularities can be accounted for by the Bethe-Salpeter equation and are removed at the 
level of the effective action by adding a contribution <5A BS to the total counterterm SX in the eight 
diagram. 

B. General nested Bethe-Salpeter equations 



The equations of motion can always be written as 

S_i(P) = \l A^(P,Q)5D(Q) + ± I A^{P,Q)a (p {Q) + E' r (P) 



and 



A H> (P,Q)M(Q)E_ 1 (Q) 



Q 



" A iptp (K,Q)6D(Q) + l I A w (K,Q)cr v (Q) + W r (K) 



Q 



+ f A H ,{K,Q)M{Q)^- 1 {Q) 
Jo 



(44) 



(45) 
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where we have used the basic kernels introduced in section III Bl The first step is then to 
create a new kernel by iterating via a Bethe-Salpeter equation. It is the generalization 
of the kernel that we introduced in the one-loop approximation: 

T^(P,K) = A^(P,K)- f T H ,(P,Q)M(Q)A H ,(Q,K), 

= A^(P,K)- f A Hl (P,Q)M(Q)T^(Q,K). (46) 
Jq 

We can then write the equation of motion for the fermion as 

S-i(P) = \ f A^(P,Q)6D(Q) + ± I A^(P,Q)a v (Q) + K(P) (47) 

J Q J Q 

r^(P,Q)M(Q)£_i(Q)- / / T^(P,K)M(K)A^(K,Q)M(Q)^ 1 (Q). 

Jo Jk 



'Q JQ 

Using the equation of motion (J4l| we obtain 

S-x(P) = \\ Kp tp (P,Q)SD(Q) + ^ f A^{P 1 Q)a ip {Q) + ^{P) 

- \ f f r^(P,Q)M(Q)A^(Q,K)5D(K) 
z Jq Jk 

-Iff r^(P,Q)M(Q)A^(Q,K)5D(K) . (48) 
1 Jq Jk 

We now plug £_i into the equation of motion for the scalar ()45|) : 

U (K) = \ f A (ptp (K,Q)SD(Q) + l f A W (K, Q)a tp (Q) + I%{K) 
1 Jq 1 Jq 

+ \ f f A^(K,Q)M(Q)A^(Q,P)5D(P) 
1 Jq Jp 

i 

2 
1 
2 




I f f f A^(K,Q)M(Q)T^(Q,K)M(K)A^(K,P)SD(P) 
^ Jq Jk Jp 

A ip 4K,Q)M(Q)T^(Q,K)M(K)A^(K,P)a ip (P) . (49) 

'< jp 

By introducing the kernel 

A W (K,P) = A (fip (K,P) + f A (pf (K,Q)M(Q)A^(Q,P) 

Jq 

- f f A^(K : Q)M(Q)T^(Q : R)M(R)A^(R } P) , (50) 
Jq Jr 
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the equation of motion for the scalar propagator becomes 



2 



1 



~A w (K,Q)5D(Q) + l [ A w (K,Q)(Q)a v (Q) + W;(K). (51) 



This equation can be renormalized by means of a Bethe-Salpeter equation (|12|) and leads 
to the equation ([43)1 . H"(K) is a finite function, decreasing as 1/K 2 which ensures the 
finiteness of equation (J4"3"|. 

C. Renormalization algorithm 

The only thing that we need is a rule to construct a finite T w : 

• Compute the kernel = —5 2 &/5S 2 . 

• Construct the kernel via the first finite Bethe-Salpeter equation: 



• Compute the kernels = -25 2 $/5S5D, A^, = Al and A w = 48 2 §/5D 2 and dress 
the kernel into A w by the following construction: 




(52) 



A W (K,P) = A W (K,P) 





A^{K, Q)M(Q)T n (Q, R)M(R)A^{R, P) . 



(53) 



• Construct T w via the second finite Bethe-Salpeter equation: 



r \ptpyPf a^") r^^?(-^*5 -^*) A^p^p(yP^ a") A^pip^Pxj a"^) 




This is equation is finite since the properties of the kernel A W (P, A') revealed in the one-loop 
case, still hold at higher number of loops as shown in Appendix [Bj 
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As an example we give the divergent structures generated by the inclusion of the fourth 
diagram in Fig. [TJ The contributions to A^, and the corresponding are depicted in 
Fig. 1221 T?/>v nas t° be dressed into A w thanks to the diagrams in A w and A^. These 



A 



FIG. 22: Kernel and ladder resummation T 
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diagrams are depicted in Fig. 031 together with the corresponding A w . T w contains all the 



> v 



\ / 



/ — >- 



FIG. 23: Building A w out of T^. 



ladders build from the rung A 



V. CONCLUSION 



We have shown that ^-derivable approximations are renormalizable in a model where 
scalar degrees of freedom are coupled to fermions via a Yukawa interaction. If the Yukawa 
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coupling renormalization can be understood by drawing boxes on resummed diagrams, it 
is not so for the scalar coupling renormalization which amounts to renormalizing a set of 
nested Bethe-Salpeter equations for the scalar and fermionic four-point functions. We give 
a systematic rule to construct and renormalize such nested equations at any loop order in 
the ^-derivable loop expansion. 

This analysis prepares the ground for QED since the power counting rules are strictly the 
same. Thus one is able to analyse all the UV singularities which may appear in ^-derivable 
approximations for Abelian gauge theories. In particular, there are singularities involving 
four photons. In perturbation theory, the latter are constrained to vanish by the Ward- 
Takahashi identity. However in any finite loop approximation of the 2PI effective action, 
these singularities are present since the different channels necessary for the cancelation to 
occur are not resummed at the same time. This is for example clear with the one-loop 
topology of Fig. ^2 which arises either from the second two-loop diagram in Fig. ^ (after 
cutting one photon line, iterating the equations of motion once and cutting a second photon 
line, see Fig. EJ) or from the fourth three-loop diagram in the same figure (by directly cutting 
two photon lines). The way to remove consistently all these spurious singularities will be 
presented in a subsequent work 
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APPENDIX A: GENERALIZED SPIN INDICES 

Let us consider a specific example, namely that of the Bethe-Salpeter equation (fTU|) . 
Writing all the indices, it reads 

r% m " s \p,K) = a^ 6 \p,k)- I A^ ) ' ( ^'Hp,g)^ v (g)r^^W(g,K)^(Q). 

JQ 

(Al) 

This ordering corresponds to the standard way of writing the diagram by following the 
fermion lines. We call it cycle since it first goes up in the diagram and then goes down, 
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following the fermion orientation backwards. There is a more appropriate way to organize 
indices in the case of Bethe-Salpeter equations. We call this ordering ladder. It reads 

r { tf' m {p,K) = a% p) ' m (p,k)- [ Ay' (a '^Hp,g)^ y (g)^ /3 ,(Q)ri 7 ;' ) ' (75) (Q,^). 

JQ 

(A2) 

Now if we define the four index object M( a / J g/) ) (y < 5/)(Q) = S a iy(Q)Ss'p>(Q) and if we interpret 
the pair of indices between brackets as generalized indices, we can rewrite the previous 
equation in terms of a matrix product in this generalized index space. We thus have 

T^(P,K) = A^(P,K)- [ A^(P,Q)M(Q)T^(Q,K). (A3) 

JQ 

This matrix notation is perfectly suited for the ladder structure of the Bethe-Salpeter equa- 
tion. For instance, iterations of the latter amount to matrix multiplication on the right. 
This idea can be extended to other situations. For example, if we consider the one fermion 
loop contributing to the scalar four-point function (see Fig. I1(J|) . we have to consider the 
trace 

/ = / S a ,(Q)A^\Q,P)S 5 ,(Q)Af;\Q,K) 
JQ 

= I A^\Q,P)M {S , UM (Q)A^\Q J K). (A4) 

JQ 

By interpreting A^(Q,P) as a vector in the generalized index space and by defining the 
transposed vector by (A^)^)(P,Q) = A£f(Q,P), we can write J in a compact form as 

/= / A^{P,Q)M(Q)A^(Q,K), (A5) 

JQ 

without using the trace symbol "tr". By definition of the kernels (see Eq. (|8J)), we have 
A Vv = A^ so that 

/ = / A^(P,Q)M(Q)Ki xp (Q,K), (A6) 

JQ 

APPENDIX B: UV PROPERTIES OF THE KERNEL A 

We describe here the general properties of the kernel A w at any loop order. A w is 
constructed from A^ in two steps. First one iterates A^, through the Bethe-Salpeter 
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equation and generates T^. Next one combines together with A w , and A^ 
according to Eq. (jo^j) . 

If the approximation belongs to the class of those described in section II VI the BPH 
technique allows to renormalize A^. Equation (J52j) is then automatically finite since, in the 
Bethe-Salpeter ladders that it generates, one cannot find subdiagrams with a non-negative 
superficial degree of divergence 5. It follows that is finite. This kernel T^, when 
combined with A^ and A^ v in order to build A w , generates divergent structures, since the 
overall degree of divergence is 5 = 0. However there are again no divergent subdiagrams. 
A w contains only an overall divergence and any combination such as A W (P, Q) — A W (K, Q) 
is finite. 

Finally, at fixed P and K, the difference A W (P, Q) — A VV> (K,Q) decreases like 1/Q at 
large Q, that is faster than each of the terms of the difference, which go like InQ. The 
reason for this is that the leading logarithmic behavior of A W (P, Q) at large Q and fixed 
P is independent of P and thus vanishes in the difference. To understand this result one 
applies Weinberg's theorem The leading logarithmic behaviour at large Q and fixed 
P is related to the subgraphs that one can draw on A w which contain the external legs 
attached to Q and whose degree of divergence is 5 = 0. Since the two legs attached to Q 
are scalar lines and since 5 = 0, the subgraph must have two other scalar external legs. The 
problem consists then in enumerating all the possible subgraphs with four scalar legs (two 
of them being attached to Q) that one can draw on A w . Since A w is 2PI with respect 
to scalar lines, the only possible subgraph is the diagram itself. The leading logarithmic 
behavior of this graph at large Q and fixed P is obtained by considering the regime where 
the momenta of the lines of the graph go to infinity simultaneously. This washes out the 
presence of P in any of the lines of the diagram and thus the leading asymptotic behavior 
is independent of P. 

APPENDIX C: THERMAL SEPARATION 

We explain here how to separate the explicit thermal dependence in the case of the 
one-loop example considered in section ITTT1 
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1. Fermionic self-energy 



We consider first the Matsubara sum appearing in the equation of motion for the fermionic 
propagator: 

/ D(Q)S{Q + P) = -Y, / D(Q)S(Q + P), (CI) 

JQ P „ Jq 



with J q = J (^d-i • m our notations Q = (iu n ,q) and P = (iu>,p), iu n and ioo being 
scalar and fermionic Matsubara frequencies respectively. We now introduce the spectral 
representation for each of the propagators: 

M r o,p + q) 



D( Q )= [ ^a), 5(g + p)= / 
J qo qo - ™ n J ro 



r — iuj — iuJ n 



with f = f + °° ^. Thus 



9o 

(T) 



£(Q)S(Q + P) = / / / P¥ ,(g , q ^(r ,p + q)l V — 

iq Jqo Jr V V 



- iu n )(r -iu - iu n ) ' 
(C2) 

We now perform the Matsubara sum: 

„ (Qo - iv n )(r - ioj - ioJ n ) r - q - iu r Q - q Q - iu ' 
The explicit thermal dependence is extracted by means of 

^90 = -0(-?o) + 4%) N M > n -r = 0(r ) - e(r )n| ro | . (C4) 



Plugging this decomposition into Eqs. (jC3|) and ()C2|) . the terms with 6 1 functions can be 
rewritten as 

D{Q)S(Q + P), (C5) 



and the terms explicitly dependent on temperature as 



/ / e{ro)pip{ro,p + q)n\ ro \D(ro-iu,q)+ / / e(q )p v {q , q)N\ qo \S{iu> + g„, P + q) • (C6) 

JqJro J q g 

After changing variables (r ,p + q) — > (— qo, — q) and using p^,(— ?o, — q) = Pip(qo,<l) and 
D(-q , -q) = D(5o, q), we obtain 

-// e(g )p^(go,q)^| go | £, (^ + 9o,P+q)+ / / e(go)p^(go,q)A r ko|' 5 (^ + 9o,P + q) • (C7) 

Jq Jan Jq Jon 
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Introducing the notations 

0>(?o, q) = e(?o)% |M?o, q) , ^(?o, q) = e^o^M/V^o, q) , (C8) 

and designating by Q = (go, q) Minkowskian momentum, we can finally write 

/ D(Q)S(Q + P) = [ D(Q)S(Q + P) 
Jq Jq 

+ / S(P + Q)a v {Q) - [ a f (Q)D(Q + P). (C9) 

2. Scalar self-energy 

The Matsubara sum involved in the scalar self-energy reads 

»(T) 



/ S t (Q)S(Q + K) = -J2 PiQWQ + K). 
Jq P n J* 



(CIO) 



After using the spectral representation, we obtain: 

' T) S\Q)S(Q + K)= I I /"^(g ,q)p^(r ,q + k)i]r 

Jrn J an J a P „ 



(g - iw n )(ro -iw - iu; n ) ' 

(Cll) 

The Matsubara frequency ia> n is now fermionic and we obtain 

\ _ ~ n qo + n rp (C12) 

f3„( q °~ iu n){r -iu - iu n ) r - g - iu ' 
Extracting the thermal dependent part from the thermal factors, we obtain for the thermal 
contribution 

- S t (r -iuj,q)a f (ro,q+k) - / o$>(q , q)S(iu + q , q + k) . (C13) 

Jro Jq Jqo J q 

Using a simple change of variables together with the parity properties of and S, one 
arrives at 

/ ( } S t (Q)S(Q + K) = [ S t (Q)S(Q + K) - 2 / a^{Q)S(Q + K) . (C14) 
Jq Jq Jq 

3. Remark on a^p 

We finish by a remark on a^(qo, q)- It is a matrix that we can decompose into rotation 
invariant components. In the massless case the decomposition reads 

<ty(?o, q) = a(?o, vho + %o, g)qV , (C15) 
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where a and b are respectively odd and even functions in go- Now if we integrate over Q, we 
obtain 

/ M<?o,q) = 7o / a(qo,q) + i l / 6(g ,g)q J = 0. (C16) 

Jq Jq Jq 
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